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Abstract 

The Stokes-Einstein relation, relating the diffusion and viscosity coefficients D and ry, is tested 
in two dimensions. An equilibrium molecular-dynamics simulation was used with a Yukawa pair 
potential. Regimes are identified where motion is diffusive and D is meaningful. The Stokes- 
Einstein relation, Dr] oc ksT, was found to be violated near the disordering transition; under 
these conditions collective particle motion exhibits dynamical heterogeneity. At slightly higher 
temperatures, however, the Stokes-Einstein relation is valid. These results may be testable in 
strongly-coupled dusty plasma experiments. 
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Two-dimensional systems in crystalline or liquid states ^ are of interest in various fields 
of physics. Monolayer particle suspensions can be formed in colloidal suspensions 3] and 
dusty plasmas 0]. Electrons on the surface of liquid helium form a 2D Wigner crystal Q]. 
Ions in a Penning trap can be confined as a single layer of a one-component plasma (OCP) ^J. 
Magnetic fiux lines in 2D high-temperature superconductors form patterns of hexagonally 
correlated vortices 0. Granular fiows of discrete macroscopic particles can be confined to 
a 2D layer 0. Here we will be concerned with 2D liquids, with no motion in the vertical 
direction. The experimental system listed above actually have finite out-of-plane motion, 
our result will be most applicable to systems such as dusty plasma where the out-of-plane 
displacement is small compared to the particle spacing. As in 3D liquids 8], 2D liquids 
near the disordering transition are especially interesting because they exhibit dynamical 
heterogeneities, where the collective motion of the more mobile particles occurs in strings . 

Our 2D system is composed of molecules or particles that interact with a Yukawa_pair 
potential. This is a soft potential that is applicable in several fields including colloids 
dusty plasmas 0, [l^ , and some polyelectrolytes [ill in biological and chemical systems. 

The Stokes-Einstein relation for 3D liquids 



apair 



D = kBT/cTiriR (1) 

is an important hydrodynamic law relating the diffusion coefficient D of a Brownian particle 
and the fiuid shear viscosity rj. It states that the diffusion of a Brownian particle in a liquid 
is proportional to the liquid temperature ksT and is inversely proportional to the viscosity 
of the liquid. This combines the Einstein relation D = ksT/^ for D and the Stokes law 
C, = CT]R for the frictional force ^ on a sphere in a fluid. Here, c is a constant that depends on 
the hydrodynamic boundary condition at the particle surface, and R is the effective radius 
of the Brownian particle. 

For 3D liquids, the validity of the Stokes-Einstein relation is well known, even at the 
molecular level. Over a wide range of .e.pe.atu.e Q. .he di<I„.on coefficient O and 
the factor kBT/r] have a flxed ratio as ksT varies. Only for deeply supercooled liquids 
does the Stokes-Einstein relation appear to fail such liquids typically exhibit a spatial 
heterogeneity [^. The Stokes-Einstein relation is so trusted in 3D that it is routinely used 
to compute D from a measurement of t]. However, the situation in 2D is very different. 

For 2D liquids, unlike 3D liquids, tests of the validity of the Stokes-Einstein relation are 
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lacking. Indeed, some earlier results suggest that the Stokes-Einstein relation should not 
be valid at all in 2D because the transport coefficients D and 11 are themselves not valid. 

n 

These earlier results came from a 2D molecular dynamics (MD) simulation [1J| and kinetic 
and mode-coupling theories 15, 2. They predict that, unlike the case of 3D, the long-time 
tails of the velocity (VACF) and stress autocorrelation functions (SACF) have a slow 1/t 
decay. A decay this slow precludes the use of Green-Kubo relations to calculate D and 
?7, and suggests that these transport coefficients are meaningless, at least in the context of 
equilibrium liquids. However, there is no universal agreement that D and 77 are meaningless 



in 2D. The shear viscosity coefficient in 2D has now been shown to 
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)e meaningful, at least for 
and an experimental 



soft potentials , in a significant number of simulations 
measurement j2l| of a 2D viscosity has been reported. The diffusion coefficient in 2D, on 
the other hand, has not been demonstrated to be meaningful over a wide range of density 
or temperature. As a test of this, motion is deemed to be diffusive, and diffusion coefficients 
are meaningful, only if the VACF decays faster than 1/t, or if the Einstein relation MSD oc t 
is obeyed. Here, MSD is the time series of mean-square displacement of particle positions. 



Simulations 2^|2^|2^, using a r~^^ pair potential and the Einstein relation, for example, 
successfully yielded a diffusion coefficient, but only near disordering transition. Identifying 
the conditions where a 2D diffusion coefficient is meaningful is our ffist goal. We do this two 
ways, using the Green-Kubo method and the Einstein relation. The expressions we used for 
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the VACF and Green-Kubo relations were Eqs. (7.2.1) and (7.2.8), respectively, in Ref. 
where we replace the coefficient 3 with 2, to reflect the number of dimensions. 

Our second goal is to test the Stokes-Einstein relation for 2D liquids. For the regimes 
where D is valid we will compute the product Drj and plot it versus temperature as an 
empirical test of the 2D Stokes-Einstein relation, 

D = kBT/c2DT^ri. (2) 

Equation Q is analogous to Eq. ((T)), except that it does not contain the particle radius R in 
the denominator due to the different dimensionality of 77 in 2D. Finding a linear scaling of 
Dt] on ksT would be an empirical verification of the validity of the Stokes-Einstein relation, 
while a deviation from linear scaling would indicate a violation. In performing this test, it 
is important that we obtained values of D and 77 independently. 

We use a 2D Yukawa system. Particles are present only in a single monolayer, and they 
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interact with a Yukawa potential energy, which for two particles of charge Q separated by 
a distance r is U{r) = {Q"^ /Aireor) exp{—r / Xo)- This potential changes gradually from a 
long-range Coulomb repulsion to a hard-sphere-like repulsion as the screening parameter 
K = a/Xf) is increased. Here, A/) is a screening length, a = (nvr)^^/^ is the 2D Wigner- 
Seitz radius and n is the areal number density of particles. A Yukawa system is 
characterized by two dimensionless parameters k and the Coulomb coupling parameter F = 
Q'^/AneoakT. We will use F"^ as a dimensionless temperature. The phonon spectrum of 
a 2D Yukawa system is characterized by a frequency Upd = ((5^/27reoma^)^/^ 26 1. The 
disordering transition for a 2D Yukawa system occurs at F = 137 for k = \4. At our value 

n 

of K = 0.56, the transition is at F ?a 145 j27|], and we denote the corresponding melting 
temperature as T^. 

We performed an MD simulation to calculate r] and D. The equations of motion for 
1024 or 4096 particles were integrated using periodic boundary conditions. A Nose-Hoover 
thermostat was applied to achieve a constant T. As a test, we ran simulations with the 
thermostat turned on or off, and verified that the thermostat does not affect rj or D, for 
any values of F. The screening parameter was k = 0.56. We recorded particle positions and 
velocities, and then we calculated velocity and stress correlation functions and the mean- 
square displacement. The system was in equilibrium, and it had no macroscopic velocity 
shear or density gradient. There was no friction term in the equation of motion. Examples of 
particle trajectories for F = 104 and 129 are shown in Fig. H Data for MSD and correlation 
functions were averaged over multiple runs with different initial conditions. Further details 

n 

of the simulation can be found in Ref. jl7l |. 

To meet our first goal, we identify regimes where the transport coefficients D and t] are 
meaningful. Previously [17[ we found that i] was meaningful, as indicated by a SACF that 
decayed faster than 1/t, for a wide range of F. Here, we report results for the diffusion 
coefficient D, where the results are more complicated. We applied two tests to determine 
whether motion is diffusive. First, if MSD oc is satisfied at times that are long (compared 
to the short-term ballistic motion, which has an exponent of 2, MSD oc t^), then the Einstein 
relation is satisfied, indicating diffusive motion. Results for this test for selected values of F 
are shown as a time series for MSD in Fig. |21 The time series has a duration that exceeds, 
by approximately a factor of 2, the ratio of the box width and sound speed. We test the last 
portion of this time series, for an interval of ~ 37ujpJ, to determine whether the exponent 
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is equal to 1. Second, if the long-time tail in the VACF decays faster than then Green- 
Kubo theory indicates that the motion is diffusive. We found that these two tests always 
agreed. Details of the VACF data will be reported elsewhere. 

We identified four different temperature regimes, and we characterized them according 
to whether the motion is diffusive or not. All four regimes have F < 145, i.e., they are 
all in the liquid state and have T > T^. In regime I (124 < F < 145) the motion is 
diffusive; this regime is near the disordering transition. Motion is also diffusive in regime 
II (88 < F < 124), and the non-ideal gas regime IV (F < 5). However, motion was not 
diffusive in regime III (5 < F < 88). That a regime exists where motion is non-diffusive is 
not surprising, because earlier studies of 2D liquids predicted they would be non-diffusive 
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due to the long-time tail in the VACF, which arises from hydrodynamic modes 
What is surprising is that we found other regimes, for our frictionless atomic system, that 
are diffusive. 

If a friction term were added to the equation of motion to model drag due to a solvent 

n 

the motion might more likely be diffusive. Lowen [2^ used a different simulation method, 
Brownian dynamics, for a 2D Yukawa system with a large friction. The equation of motion 
in that method lacks an inertial term, and is therefore suitable for a strongly dissipated 
system. They reported finding diffusive motion not only in our regimes I and II, but also for 
one data point in our regime III. Their values for D were systematically smaller than ours 
by a factor of 2.3 in regime II. 

Our two methods of computing D yield values that agree, as shown in Fig. |S1 Note that 
D increases with temperature, and the rate of increase is very high in regime I, near the 
disordering transition at T^. 

To meet our second goal of testing the Stokes-Einstein relation, we compute the product 
Dt] and examine its scaling with temperature. Results, shown in Fig.HJ reveal that Dt] varies 
almost linearly with temperature in regime II, but not regime I. We therefore conclude that 
the Stokes-Einstein relation is obeyed in regime II but not in regime I. 

We can thus summarize our two chief results. First, the diffusion coefficient is meaningful 
in both regimes I and II but not regime III. Second, the Stokes-Einstein relation is valid in 
regime II, but it is violated in regime I. 

We found the value of the coefficient in Eq. Q relating D, rj, and ksT. Fitting the data 
in Fig.|3]to a straight line passing through the origin yields C2d = 1.69±0.10. This coefficient 
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has never been reported, to our knowledge, for 2D, although its value is well-known for 3D, 
where c = 4 for stick or c = 6 for slip boundary conditions. 

Verification of the Stokes-Einstein relation in regime II indicates a coupling of diffusion 
and shear viscosity in 2D just as it does in 3D simple liquids. For a macroscopic particle 
in a fluid, particle diffusion is a result of Brownian motion, with a friction that arises from 
the viscous dissipation of the associated velocity field in the surrounding liquid. In a simple 
liquid, validity of the Stokes-Einstein relation is an indication of a coupling of individual 
particle motion to collective modes. This coupling has been observed previously in 3D for 
strongly-coupled Yukawa systems j3] and one-component plasmas [2^. 

The violation of the Stokes-Einstein relation that we observed in regime I, near the 
disordering transition, indicates a decoupling of diffusive and viscous transport. We found 
the diffusion coefficient is larger than would be expected if the Stokes-Einstein relation were 
valid. We note that in 3D a violation of the Stokes-Einstein relation is also observed in 
supercooled liquids 3,1^1, and this motivates us to suggest an explanation for our result 
in 2D. 

The decoupling in regime I (near the disordering transition) is likely related to structural 
and dynamical heterogeneities in the liquids. As shown in Fig.Q on a time scale of tens of ujpd 
the 2D liquid develops domains of localized crystalline order with localized oscillatory motion 
that are separated by stringlike structures of more mobile particles. For 3D supercooled 
liquids, it is believed that this type of heterogeneity accounts for the failure of the Stokes- 
Einstein relation Our results suggest that the same effect is responsible in 2D. 

The results we have reported for our test of the Stokes-Einstein relation in a 2D liquid 
could be tested in a future experiment. Monolayer suspensions of microsphere in a dusty 
plasma have a Yukawa interaction as in our simulation, and a recently-developed method 
has been reported for measuring the shear viscosity. There is no obstacle to measure the 
self-diffusion coefficient as well. 

We thank V. Nosenko and F. Skiff for helpful discussions. This work was supported by 
NASA, DOE, and CRDF No. RU-P2-2593-MO-04. 



[1] K. J. Strandburg, Rev. Mod. Phys. 60, 161 (1988). 

[2] C. A. Murray, W. O. Sprenger, and R. A. Wenk, Phys. Rev. B 42, 688 (1990). 

6 



[3] S. Nunomura, D. Samsonov, and J. Goree, Phys. Rev. Lett. 84, 5141 (2000). 

[4] C. C. Grimes and G. Adams, Phys. Rev. Lett. 42, 795 (1979). 

[5] T. B. Mitchell et al, Phys. Plasmas 6, 1751 (1999). 

[6] P. G. Gammel et al, Phys. Rev. Lett. 59, 2592 (1987). 

[7] H. M. Jaeger, S. R. Nagel, R. P. Behringer, Rev. Mod. Phys. 68, 1259 (1996). 

[8] S. C. Glozer, J. Non-Cryst. Solids 274, 342 (2000). 

[9] C. Reichhardt and C. J. O. Reichhardt, Phys. Rev. Lett. 90, 095504 (2003). 

[10] O. S. Vaulina and O. F. Petrov, JETP 99, 510 (2004). 

[11] A. R. Denton, Phys. Rev. E 67, 011804 (2003). 

[12] H. J. Parkhurst and Jr. and J. Jonas, J. Chem. Phys. 63, 2705 (1975). 

[13] E. Rossler, Phys. Rev. Lett. 65, 1595 (1990). 

[14] B. J. Alder and T. E. Wainwright, Phys. Rev. A 1, 18 (1970). 

[15] M. H. Ernst, E. H. Hauge, and J. M. J. Van Leeuwen, Phys. Rev. Lett. 25, 1254 (1970). 

[16] J. R. Dorfman and E. G. D. Cohen, Phys. Rev. Lett. 25, 1257 (1970). 

[17] B. Liu and J. Goree, Phys. Rev. Lett. 94, 185002 (2005). 

[18] M. Ferrario, A. Fiorino, G. Ciccotti, Physica A 240, 268 (1997). 

[19] W. G. Hoover and H. A. Posch, Phys. Rev. E 51, 273 (1995). 

[20] D. Gravina, G. Ciccotti, and B. L. Holian, Phys. Rev. E 52, 6123 (1995). 

[21] V. Nosenko and J. Goree, Phys. Rev. Lett. 93, 155004 (2004). 

[22] R. Zangi and S. A. Rice, Phys. Rev. Lett. 92, 035502 (2004). 

[23] D. N. Perera and P. Harrowell, Phys. Rev. Lett. 81, 120 (1998). 

[24] M. M. Hurley and P. Harrowell, Phys. Rev. E 52, 1694 (1995). 

[25] J. -P. Hansen and I. R. McDonald, Theory of Simple Liquids, 2nd ed. (Elsevier Academic 
Press, 1986). 

[26] G. J. Kalman et al, Phys. Rev. Lett. 92, 065001 (2004). 

[27] P. Hartmann, G. J. Kalman, Z. Donko, and K. Kutasi, Phys. Rev. E 72, 026409 (2005). 

[28] H. Lowen, J. Phys. Condens. Matter 4, 10105 (1992). 

[29] H. Ohta and S. Hamaguchi, Phys. Plasmas 7, 4506 (2000). 

[30] P. Schmidt, G. Zwicknagel, P. -G. Reinhard, and C. Toepffer, Phys. Rev. E 56, 7310 (1997). 

[31] M. Dzugutov, S. L Simdyankin, and F. H. M. Zetterling, Phys. Rev. Lett. 89, 195701 (2002). 



7 



FIG. 1: Particle trajectories from the simulation. Regime I is near the disordering transition, 
and regime II is at a slightly higher temperature. Regime I has a dynamical heterogeneity; strings 
of mobile particles flow between domains consisting of less-mobile caged-particles. In regime II 
(r = 104) the domains are smaller and there is a less distinctive difference between the caged 
trajectories and the trajectories of particles bounding the domains. Trajectories are shown for a 
time interval cOpd^t ~ 37. Only a portion of the simulation box is shown here. The full simulation 
box was 113.98a x 98.16a with 4096 particles for T = 129, and 56.99a x 49.08a with 1024 particles 
for r = 104. 



FIG. 2: Mean-square displacement (MSD). Ballistic motion at small t has MSD oc t'^, with a slope 
of 2 on this log-log plot. Motion at later times is either diffusive, MSD oc t as indicated by a slope 
of 1, or non-diffusive as in regime III. By this test, we found that motion is diffusive in regimes I, 
II and IV, but not III. A separate test, using the VACF, is consistent with this result. 



FIG. 3: Diffusion and viscosity coefficients. Diffusion coefficients (open symbols) were computed 
from the same simulation data using two methods, the Green-Kubo relation and the Einstein 
...io. Fo. . .he .ff„sio„ coefflcie., ... .sco* coefflcie. .a.a .o. U., Q 

are shown as solid circles; these were computed using a Green-Kubo relation. The disordering 
transition for our k = 0.56 occurs at F w 145, based on the simulation results of Ref. 



FIG. 4: Test of the Stokes-Einstein relation. Linear scaling Drj oc T in regime II indicates the 
validity of the Stokes-Einstein relation. However, the deviation from the linear scaling for F > 124 
indicates a violation of the Stokes-Einstein relation in regime I, which is near the disordering 
transition. Diffusive motion was not observed in regime III, so that the relation is not testable 
there. The vertical axis is normalized by nmoj'^^a^. Extrapolating the Stokes-Einstein line shown 
here into regime IV (not shown), we found rough agreement with simulation data even though the 
applicability of the Stokes-Einstein relation in non-ideal gas is not expected. 
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This figure "figure l.png" is available in "png" format from: 



http://arXiv.org/ps/cond-mat/05 1 1209vl 
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